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Abstract 

The tremendous capacity gains promised by space division multiple access (SDMA) depend critically on the 
accuracy of the transmit channel state information. In the broadcast channel, even without any network interference, 
it is known that such gains collapse due to interstream interference if the feedback is delayed or low rate. In this 
paper, we investigate SDMA in the presence of interference from many other simultaneously active transmitters 
distributed randomly over the network. In particular we consider zero-forcing beamforming in a decentralized (ad 
hoc) network where each receiver provides feedback to its respective transmitter We derive closed-form expressions 
, for the outage probability, network throughput, transmission capacity, and average achievable rate and go on to 

quantify the degradation in network performance due to residual self-interference as a function of key system 
parameters. One particular finding is that as in the classical broadcast channel, the per-user feedback rate must 
' increase linearly with the number of transmit antennas and SINR (in dB) for the full multiplexing gains to be 

preserved with limited feedback. We derive the throughput- maximizing number of streams, establishing that single- 
stream transmission is optimal in most practically relevant settings. In short, SDMA does not appear to be a prudent 
design choice for interference-limited wireless networks. 

O 

I. Introduction 

In multiuser MIMO (multiple-input, multiple-output) channels, the spatial multiplexing capability of- 
^ fered by multiple antennas can be advantageously exploited to significantly increase the achievable 
\^ throughput. In single-cell point-to-multipoint channels, the achievable sum rate scales linearly with the 
number of transmit antennas, even when the mobile users have only a single antenna. By duality, linear 
increase with the number of receiver antennas can also be achieved in multipoint-to-point channels, even 
^ with a single-antenna transmitter. Extensive research on MEMO broadcast channels over the last few years 
^ has revealed that the capacity can be boosted by transmitting to multiple users simultaneously, by means 
^ ! of Space Division Multiple Access (SDMA), rather than trying to maximize the capacity of a single- 
> ■ user link [[T|, [|2l. Nevertheless, all these promising gains critically depend on accurate channel state 
^ information (CSI), and in contrast to point-to-point channels, the quality of CSI affects the multiplexing 
^ gain of multiuser MIMO systems. As a result, a considerable amount of effort has been dedicated to 
multiuser MIMO systems operating with partial CSI at the transmitter (CSIT) in the absence of out-of- 
cell interference [|3l. 

In this paper we are interested in the capacity gains that SDMA may provide in decentralized (ad hoc) 
networks with both intra-cell due to imperfect CSIT and other user interference due to uncoordinated 
concurrent transmissions. We aim at answering whether and how aggressive use of multiple antennas 
through SDMA may increase the network throughput under a broad set of scenarios. Specifically, we 
build upon the practically relevant limited feedback model, in which each user is allowed to feed back 
5-bit quantized information on its channel direction based on a predetermined codebook known at both 
the transmitter and the receivers. Due to the high implementation complexity and sensitivity to channel 
errors of the optimal scheme (dirty paper coding [4]), linear precoding based on zero-forcing beamforming 
(ZFBF) [m is employed here. ZFBF has been shown to achieve full multiplexing gain while exhibiting 
reduced complexity [6J, [|7J. ZFBF is also used both for its asymptotic optimality in MIMO point-to- 
multipoint channels and its analytical tractability. In short, it is a logical starting point for understanding 
the effect of imperfect CSIT in interference-limited networks. 
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A. Related Work 

Several recent papers have studied MIMO ad hoc netvv'orks with Poisson distributed interferers, but the 
majority of prior work considers only point-to-point ad hoc links, in which each transmitter communicates 
with only one receiver at a time (TDM A). Several receive antenna processing techniques have been 
investigated quantifying the performance gains of antenna combining and interference cancelation [8 |- 
IfTTTl . Open-loop spatial multiplexing with linear receivers is studied in [|T2]| . while [|T3l studies multi-mode 
precoding with receive interference cancelation. Spatial multiplexing with limited transmit CSIT have 
been studied for asymptotically large number of antennas in [14J. Multiuser transmission, in which each 
transmitter sends different messages to multiple users, is analyzed recently for scalar broadcast channels 
using superposition coding in [fTSl . The performance of multiuser MIMO communication in a Poisson 
field of interferers is first studied in [16] considering non- linear and linear precoding with perfect CSIT. 
Therein the network-wide capacity is shown to increase linearly and sublinearly with the number of 
antennas when DPC and ZFBF are employed, respectively, for single-antenna receivers. It is also shown 
that it is often throughput maximizing to send fewer streams than the number of transmit antennas. Here, 
we investigate whether these results derived under the idealized perfect CSI assumption still hold in a 
more practical and realistic scenario. Single-user beamforming is compared to SDMA with perfect CSIT 
in ifm for two-tier networks with spatial randomness. A key finding is that single-stream transmission at 
each tier provides significantly superior coverage and spatial reuse relative to multiuser transmission. 

B. Main Contributions 

In this paper, in contrast to all prior work on multiuser MIMO ad hoc networks, we consider SDMA 
ad hoc links in which only local, partial CSI is known at the transmitter. We investigate the zero-forcing 
transmission technique with quantized CSIT in the context of decentralized interference-limited networks 
using a random access medium access control protocol. The spatial distribution of the nodes follows 
a homogeneous Poisson point process (PPP), and each transmitter serves multiple receivers (point-to- 
multipoint), each located at a certain distance away from it. 

First, we investigate the performance of limited feedback SDMA and derive novel closed-form ex- 
pressions for the outage probability, network aggregate throughput, transmission capacity [18J, and mean 
user rate using stochastic geometric tools. These expressions enable us to quantify the capacity gains that 
zero-forcing precoding can provide in decentralized networks and the network-wide performance loss due 
to residual (uncanceled) multiuser interference caused by quantized channel information. Key findings are 
that there is a node density value that maximizes the network throughput and that for practically relevant 
values of feedback bits, if the number of antennas/streams is increased, the outage constraint for non zero 
node density is prohibitively high. 

Second, we evaluate the performance degradation as a function of the feedback rate, the number of 
antennas, and the outage constraints for fixed feedback levels. A first result is that, similar to the single-cell 
case, the throughput of feedback-based zero-forcing is self-interference-limited, i.e. the average throughput 
is bounded if the feedback load is kept fixed, even if the transmit power is taken to infinity. We then 
provide the scaling of feedback bits in order to guarantee bounded performance offset. A key finding 
is that per-user feedback load B must be increased almost linearly with the number of antennas and 
logarithmically with the target signal-to-interference-plus-noise ratio (SINR), i.e. similar to the single-cell 
case analyzed by Jindal [19] if the transmit power is replaced with the target SINR constraint. 

Finally, we derive the optimal number of streams in order to maximize throughput and transmission 
capacity, enabling us to propose transmission schemes that dynamically adjust the number of streams to 
the network operating values. Our results establish that the optimality of single-stream transmission in 
most practically relevant settings. The main takeaway of this paper is that limited feedback SDMA may 
not be a wise use of an antenna array in ad hoc networks, and in general in network settings where 
interferers can be close by as in heterogeneous networks including femtocells, picocells, relays, and WiFi 
hotspots. 
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II. System Model and Preliminaries 

The network model consists of transmitters arranged according to a homogeneous Poisson point process 
(PPP) $ of intensity A in M^. Each transmitter has M antennas and communicates with a set of intended 
single-antenna receivers /C with cardinality \}C\ = K < M. Every transmitter sends K streams destined 
to a different user each, forming thus a i^-user broadcast cluster. Users are distributed according to some 
independent stationary point process. Transmissions are uncoordinated (random access MAC protocol) 
and the signal is attenuated according to the standard power law, i.e. the received power decays with the 
distance d as rf^" for path loss exponent a > 2. For the fading model (random channel component), we 
assume that all point-to-point channels experience i.i.d. Rayleigh block fading with unit mean. 

We investigate multiuser transmission, in which each transmitter sends a different message to its associ- 
ated receivers. Due to the stationarity and the rotational invariance of the PPP (Slivnyak's Theorem [|20l ). it 
is sufficient to analyze the performance of a typical one-to-many communication channel, as the statistics 
of the signal reception remain the same for the disjoint set of intended receivers in the broadcast cluster. 
Denote the typical transmitter by Tq located at the origin communicating with its k-ih typical user, located 
at distance d^, for G /C and denoted as R^q\ 

The received signal y^k at typical receiver R^q \ k G JC assuming frequency-flat channels is given by 

Vok = \fp dj^^ hofcXo + y/p -Dj ' hifcXj + riofc (1) 

je'i'(A)/To 

where p = ^, P is the transmit power, Di E M? is the distance to the i-th transmitter, and riofc is the 
complex additive Gaussian noise with variance a^. The vector channel from the i-th transmitter to R*^^ 
is denoted by hik E C^^^^ and Xj is the M x 1 normalized transmit signal vector of the i-th transmitter. 
The vectors h are assumed to have i.i.d. CJ\f{0, 1) entries, independent across transmitters and of the 
random distances Di. The index '0' is dropped for notation simplification, as all the subsequent analysis 
is performed on a typical broadcast cluster. 



A. Finite Rate Feedback Model 

We assume that each receiver has perfect knowledge of the channel to its corresponding transmitter. 
In each cluster, the transmitter acquires partial CSIT that only captures the spatial direction information 
of the channel, referred to as channel direction information. For that, a quantization codebook Vk = 
{vfei, Vfc2, • • • , VfcAf} containing N = 2^ unit norm vectors {^^ki}iLi ^ C*^ is employed, assumed to be 

(k) 

known to both Tq and receiver i?Q . At each feedback reporting slot, each receiver k quantizes its channel 
realization to the closest codeword with respect to the chordal distance [|2T|| . [|22ll . 

hfe = arg max |hfeVfci|^ = arg max cos^(Z(hfc, v^i)), 

where = h^/ ||hfc|| corresponds to the channel direction. Each user sends the corresponding quantization 
index back to the transmitter using B = \\0g2 N] bits through an error and delay-free feedback channelQ. 
The optimal vector quantization strategy in multiuser downlink channels is not known in general, even 
in single-cell systems, and is out of the scope of our work. We resort hence to a vector quantization 
scheme following the quantization cell approximation (QCA) [i22l . [|23l . It has been shown in ll23l . Il24l 
that QCA can facilitate the analysis and provide a very accurate performance approximation, with only 
small difference from random vector quantization. 

'The error-free assumption can be well approximated using sufficiently powerful error-correcting codes over the feedback link, whereas 
the zero-delay assumption may be valid when the processing and feedback delays are small relative to the channel coherence time. 
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B. Zero-forcing Beamforming 

In this paper, we focus on linear precoding (downlink beamforming), in which the transmit symbol 
vector X is a linear function x = ^^g^ w^s^, where Sk is the data symbol intended for the A;-th receiver 
and Wjt G C^^^^ is the unit-norm beamforming vector for user k. Specifically, users are served via 
zero-forcing for which the beamforming vectors are chosen as h^Wj = 0, ^ i, k E IC. 

Let H(/C) = hj , . . . denote the concatenation of the quantized channel vectors upon which 

zero-forcing is performed. The beamforming vectors are given by the Moore-Penrose pseudoinverse 

W(/C) = n{icy = H(/C)"(H(/C)H(/C)^)-^ (2) 

with Wfe obtained by normalizing the k-th column of W(/C). The receive SINR at the k-th typical user, 
treating interference as noise and using equal power allocation for each of the data streams, can be 
expressed as 

........ 



SINRfc 



L + L + o- 



2 



with Ip= Pl|h.fcW,f D-'^ (3) and = 5Z P\^k^^\^ ^^'^ (4) 

iG*(A)/To je/CjVfc 

where Ip is the aggregate inter-cluster interference from the Poisson field of interferers $/To and Ig is the 
intra-cluster self-interference due to the fact that zero- forcing vectors are calculated based on quantized 
CSIT. 



C. Performance Metrics 

Outage probability. A primary performance measure is the outage probability, which is defined as the 
probability that the received SINR falls below a target SINR (3, i.e. 

J^(/3, a) = P (SINR < /3) . (5) 

It can be thought of equivalently as the probability of no coverage of a user, and is evidently a continuous 
increasing function of the intensity aH In multiuser communication, different SINR statistics may be seen 
on different users (streams), resulting in a per-user outage probability J'^''\/3k, a), Wk E /C, with being 
the target SINR on stream k. 

Network throughput. The network throughput is defined as the product of the unconditioned success 
probability and the sum rate per unit area assuming that capacity-achieving codes are used. When K 
independent data streams are sent on each broadcast cluster the throughput is given by 

r = A 5^ P (SINRfc > /3fc) log2(l + /3fc). (6) 

Note that the success probability P (SINR^ > /3k) is itself a monotonically decreasing function of A. As 
the success probability is not constrained to a minimum value, the throughput-maximizing density may 
be obtained at the expense of very high outage levels. 

Multi-stream transmission capacity. Generalizing [|T8l for the case where K streams are sent by each 
source node, we define the multi-stream transmission capacity (TC) as the maximum number of concurrent 
multi-stream transmissions A^ per unit area allowed subject to an outage constraint e, i.e. 

C = KX,{l-e), (7) 

^In an SDMA setting, we can alternatively define = P(X(xo;yo) < r), where X(xo;yo) is the mutual information between xq 

and yo, and r is a certain target information rate. However, a decomposed per-stream/user outage constraint is more meaningful in SDMA 
ad hoc networks, in which each stream contains a different message. 
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where A, = sup {A : P (SINR^ < Pk) < Cfc, VA; G /C} (8) 

defines the maximum contention density for a per-stream outage constraint G (0,1). This outage- 
based metric quantifies how efficiently the network utilizes space as resource under a maximum outage 
constraint, as opposed to the network throughput that may result in high outage events. In other words, it 
calculates the maximum density of transmissions per unit area so that all K users in the broadcast cluster 
do not exceed a desired outage level e. 

Average ergodic rate. Finally, we define the average data rate (in nats/Hz) achievable by a typical user 
assuming Shannon capacity achieving modulation and coding for the instantaneous SINR to be 

7^(A,a) =E{log(l + SINR)}. (9) 

In contrast to the two aforementioned metrics, this average capacity measure presumes dynamic rate 
adaptation to the instantaneous SINR. 

In the remainder, for the sake of exposition simplicity we assume that all streams have identical stream 
outage constraint e and SINR target (5, i.e. = e and (5k = V/c G /C. 

III. ZFBF Performance Analysis 

In this section, we derive new closed-forms expressions for the network throughput, transmission 
capacity, and average achievable rate of zero-forcing precoding. 

A. Outage Probability 

Theorem 1: The outage probability for the k-th typical user in a wireless ad hoc network using 
multiuser zero-forcing with quantized CSIT is given by 

where Ck = hRt, 5 = 2'^^, and Xk = ^ V i]B { m + -,K - m - - ], 

" ^0 vW V « «/ 

with B{a,b) = /p^t'*~^(l — t^^^dt = ^p"j^|^^^ being the Beta function and T{x) = J^^ t^~^e~^dt the 
Gamma function. 

Proof: See Appendix |Al 

As expected, the outage probability is a decreasing function with the feedback bit rate B and an 
increasing function of K since > 0. Fig. \T\ shows the outage probability vs. the node intensity for 
different values of feedback load and antennas. The single- antenna, single-stream (SISO) outage probability 
is also plotted for comparison. Numerical evaluations of Theorem [T] confirm the analysis and show that 
increasing the number of antennas/streams results in unacceptably high outage probability values even in 
sparse networks (A — t- 0). 

B. Network Throughput 

Based on the above derived outage probability and the throughput definition (cf.®), when K indepen- 
dent data streams are sent on each broadcast cluster, the total number of successful bits/s/Hz/unit area 
(throughput) is given by 

2 

/.2/a. 2/- I / \T Cmax 
^ ^~>^kC,^ p-O- Ck/P (a) p-A-LxCmaxp p 

r = \Yl (1 + p^s)^-i l°g2(l + > KX log2(l + /3). (10) 

where (a) results by setting /3k = /3, Vfc G /C, and Cmax = /^"^max with (imax = maxdfc. The approximation is 

fce/c 

derived for the sake of exposition simplicity and can be seen as a lower bound on the network throughput 
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with equal target SINR for all users. In the sequel, for exposition convenience and unless otherwise stated, 
we consider that K = M users are served on each cluster. 

Remark 1: By taking the derivative of T with respect to A keeping all other parameters fixed, we see 
that throughput decreases with the node density, if 

A > (XM/3^rfLx)"' = A*. (11) 

The optimal intensity A* does not depend on the number of feedback bits, which only affects the 
amount of inter-cluster interference. However, although throughput can be maximized for A > A*, 
the number of outage events can be arbitrarily high due to intra-cluster interference. For A = A*, the 
success probability becomes ^ (i^^^'jTi-i > which means that multiuser zero-forcing beamforming with 
finite rate feedback decreases the success probability by a factor of (1 + (36)'^^^^ as compared to the 
case of point-to-point ad hoc communications. Furthermore, A* is a regularly varying function of M 

with index —2/a and for large M, the optimal network density scales as A* = 0(M^~), namely 

A* / 2 \-i 
lim 2- = ( vr/Jarf^ r(l — 2/a;) ) . As expected, the optimal contention density decreases when 

M^oo ^ ' 

M increases since increasing the number of streams sent boost the inter-cluster interference /p. The 
optimal density also decreases for /3 or rfmax increasing as the reliability requirements on the per-user 
performance become higher and harder to satisfy. 

Remark 2: In terms of feedback rate, the network throughput can be shown to be a monotonically 
increasing function with B. Focusing now on the transmit antenna configuration, after some algebraic 
manipulations, we can show that throughput is maximized if 

M* = max([£J ,1) 

where i is the nontrivial solution for M of 



(log 2) ■ + + \mc-:-^ = 1. 



M-1 ' ^ ' 1 + 135 p ^"^^"9M 

Although the partial derivative of Xm can be expressed as sum of beta and digamma functions ^^0(2;) 
since ^^j^^^^ = B{y, x){tpo{x) — ^o{x + y)), a closed-form expression for £ is hard to obtain. 
An analytical expression for M* can be found by applying the large M approximation Im ~ 7rr(l — 
2/a)M2/°. In that case M* = x", where x is the solution of the polynomial equation csx^" + C2X°^^ + 
(ci - C3 - l)x" - C2X^ + 1 = 0, where ci = ^ log 2, C2 = A7rr(l - 2/a)Cm4x, and C3 = From 
Abel's impossibility theorem fIS], a formula solution only exists for a < 4, while for a = 3 the solution 
can be expressed using Kampe de Feriet functions. 

In Fig.[2lwe plot the network throughput vs. the intensity A. We observe that throughput is a decreasing 
function of the number of antennas and that the performance degradation from imperfect CSIT is more 
pronounced for M increasing. The SISO and the perfect CSTbased zero-forcing are also plotted for 
comparison. We also see that multi-stream transmission is slightly superior in sparse networks (low A), 
but is generally outperformed by SISO. 

C. Multi-stream Transmission Capacity 

We turn now our attention to the maximum achievable throughput under bounded outage levels. 
Theorem 2: The maximum multi-stream transmission capacity of limited feedback zero-forcing precod- 
ing in random access ad hoc networks is given by 

c = ('og - ^ - Ml + m"-') . (12) 

Proof: The result follows by finding the inverse of the expression J^m{P) a) = P I min SINR^ < /3 ) 



e with respect to A, i.e. J^^^(/3,a), and substituting it in ([7]). 
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The second term in (fT2l) captures the effect of background noise on multi- stream TC, whereas the 
third term corresponds to the capacity degradation from quantized CSIT. For large M with K = M, 
the first term scales as 6(Af^^~), whereas both second and third terms scale as 0(M^~~) (for fixed 
feedback quality). This implies that the detrimental effect of residual interference from quantized channel 
information is orderwise dominant, becoming the transmission capacity limiting factor. 

Interestingly, in contrast to the cases of point-to-point and multiuser ad hoc communication with perfect 
CSI, it is not guaranteed that non-zero transmission capacity can be achieved for any feedback rate due 
to the self-interference that cannot be completely eliminated with quantized CSIT. After some algebra, 
we can show that the amount of feedback resolution Bj defining the multi-stream TC feasibility region, 
i.e. the region for which positive maximum contention density exists, is 



Bf> 



(M-l)log2 



1 

A'-l — 1 



(13) 



provided that (1 — e)e'^^^"i='''/^ < 1. The latter condition is more general and applies even to systems with 
perfect CSI as it guarantees that a non negative A exists for certain SNR = p/a^ and outage constraints. 

In Figs. [3] and m we plot the transmission capacity vs. the outage constraint and the number of antennas, 
respectively. We observe that positive transmission capacity is achieved for significantly high outage 
e, while increasing the number of streams further deteriorates the performance. In practically relevant 
scenarios where the outage constraint is kept low, single-stream transmission (TDMA) is optimal, even 
in the high resolution regime (B — )■ oo). Furthermore, we see that multi-stream transmission is beneficial 
for low number of antennas/streams and for relatively large number of feedback bits. 



D. Average Achievable Rate 

In this section, we derive the average data rate achievable by a typical receiver assuming Shannon 
capacity achieving modulation and coding for the instantaneous SINR. Note that this is not the maximum 
achievable Shannon capacity as each transmitter-receiver link treats interference as noise. 

Theorem 3: The ergodic rate in nat/s/Hz of a typical receiver in a broadcast cluster where multipacket 
transmission is employed using finite-rate feedback zero-forcing is given by 

2. 

roo Cix -C2xa 

= E{log(l + SINR)} = y^ _____dx. (14) 

where Ci = a'^R'^/p and C2 = XImR\- 
Proof: See Appendix |Bl 

For general values of « > 2, the computation of the average user rate IZk involves numerical integra- 
tion. In the interference-limited regime (cr^ — t- 0), pseudo-closed-form expressions involving generalized 
hypergeometric functions can be found, but these provide little insight on how different system operating 
parameters affect the average user rate. Therefore, we consider concise performance bounds. 

As the main focus of the paper is to investigate the relationship among feedback bit rate B and both 
inter-cluster and intra-cluster interference, we first provide the following result that shows that the average 
achievable rate with limited feedback of fixed quality converges to a finite ceiling as SNR — )• 00. 

Theorem 4: The average user rate of imperfect CSIT-based zero-forcing is upper bounded by 

< ^ + i^M-i - m) + log + dr ( ppy) V 5(m - 1) + ^ I (15) 

where = J2^=i 7 '■^ ^'^^ harmonic number and ^(■) is the digamma function. 
Proof: See Appendix O 
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The above result implies that at high SNR (P — t- oo), the user rate is bounded and the system becomes 
interference-limited no matter how many feedback bits are reported back to the transmitter. The upper 
bound in the above theorem is quite loose in general, however it was derived for demonstrating the quasi 
linear dependence of the average rate and the feedback load B. We derive now a tighter upper bound 
by applying integral inequalities directly to (fT4] ) as a means to find closed-form expression for the mean 
ergodic rate. 

Lemma 1: The average achievable user rate is upper bounded by 



7^fc < min / e-"^^i^e-"i^^^"dx ■ {A{vi))~ (16) 

1 1 a 1 

< min {U1U2C1) "i"2 (r(l + a/2))"i"2 (^^1)2(72) '^'^'^''■i ■ [A{vi)Yi , (17) 



Ul,t)l,U2,t'2 



with 1 < Ml, M2, f 1, f 2 < 00, — + — = 1, — + — = 1, and 

A{v,) = {6- 1)— 5--i?(M.i - 1, 1 - (M - l)v,) + ^m,vu2-{M-l)v,,l/5) ^^^^ 

Vi[M — 1)0 

where 2Fi{a,b,c, z) denotes the 2F1 Gauss hypergeometric function, and B{x,y) the Beta function. 
Proof: See Appendix 

For ui = U2 = vi = V2 = 2, which provides the tightest upper bound in most cases, significant 
simplification is possible for (fTSl ) since ^(2) is given in terms of log(5) and a polynomial expression in 
5. In that case, we can easily show that the upper bound on ergodic capacity scales like 0(A~"/^). The 
validity of the above bounds is verified in Fig. |5l where the ergodic capacity given by (fT6l) (cf. upper 
bound 1) and (flTI) (cf. upper bound 2) is compared with the exact average user rate for different number 
of feedback bits. We observe that the tightness of our bounds is improved when the number of antennas 
and feedback bits is increased. 

In the no noise case, the above lemma results in 



M 



;i + sx) 



l-M 



1 + X 




(19) 



which scales as 6(A^^) and is inversely proportional to the distance between transmitter and the k-th 
user. For A and B that do not depend on M and for large M, by calculating the integral, the upper bound 
in (fT9l ) is shown to be inversely proportional to M. This implies that in order to have per-user mean rate 
that does not scale with the number of antennas, A ~ M^^. Finally, in the interference-limited regime 
and for a = 4, (fT4)) admits a closed-form expression in terms of Meijer-G functions and trigonometric 
integrals. The no noise upper bound (fT9] ) approximates very well the exact ergodic rate even for moderate 
SNR values for increasing number of antennas. 

Returning now to the general case, the following easily computable result can be obtained using different 
bounding techniques. 

Lemma 2: The average user rate with finite rate-based zero-forcing satisfies 

< < (20) 



with 



7^- = -2cos(C2)Ci(C2) + ^^i^Gi^(f 



-1 



2' 



+ sin(C2)(7r-2Si(C2)) (21) 



^^'^ = p^:^y^[2cos(C2)Ci(C2)-2cos(C72)Ci((72) (22) 
- vr sin(C2) + tt sin(C'2) + 2 sin(C2)Si(C2) - 2 sin(C'2)Si(C72)] (23) 
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where G'^j^^z is the Meijer-G function, Si(x) = ^ ^dt is the sine integral, Ci{x) = 

— ^^dt is the cosine integral, and C2 = p^jj^- 

Proof: The bounds are obtained by applying Bernoulli's inequality [|26l for the integrand and evalu- 
ating the resulting integrals. Specifically, for the lower bound, we apply the inequality (1 + xY > 1 + rx, 
for X > —1, r < or r > 1, to the function g(x) = (1 + 5x)^^^^. For the upper bound, we use that 
(1 + < (1 + rx)~^ for r > 0,x > -1. ■ 
In Fig. [6l the achievable mean user rate is compared with (|23T ) (cf. upper bound) and (|2TI) (cf. lower 
bound) vs. SNR. We observe that the bounds are very tight at low M for all SNR range, while the lower 
bound becomes loose at high SNR when the number of antennas is increased. 



IV. Effect of Limited Feedback 

In this section, we analyze the effect of feedback quality on the network performance and provide 
design guidelines for the system operating points based on our analytical framework. In particular, we 
show at which rate feedback has to scale to maintain a certain bounded network capacity gap. We also 
derive the optimal number of streams/users to be employed in order to maximize the network throughput 
and the multi-stream transmission capacity. 



A. Performance Degradation due to Finite Rate Feedback 

We first provide the feedback bit scaling that guarantees constant (bounded) performance loss between 
the performance of zero-forcing with perfect CSI and that with partial CSIT. 

1 ) Transmission Capacity: The transmission capacity gap AC is defined as the difference between the 
transmission capacity achieved by perfect CSIT-based and that of limited feedback-based zero-forcing, 
i.e. AC = (Ccsi —C), where Ccsi is the multi-stream transmission capacity given by (fT2l) for i? — )► 00 
(perfect CSI). Thus, the performance degradation is given by 

In order to maintain a transmission capacity offset AC = lege, after some algebraic manipulations, we 
have that the number of feedback bits per user satisfies 

Bac > (M - 1) log2 /3 - (M - 1) loga I c^i-^-Dd-^) - 1 1 bits/user, 

K{K-1){1-,) 

with 1 < c < (1 + f3) ^KiilS^ for a non trivial result. Therefore, to guarantee a constant performance 
offset in terms of transmission capacity, the number of feedback bits per user must be increased at least 
linearly with the number of antennas/streams and approximately logarithmically with the target SINR 
constraint. This is basically the same scaling behavior as in |fT9l for fixed average rate offset if the 
target SINR constraint is interchanged with the transmit power. Interestingly, TC for ad hoc networks 
appears to capture the performance degradation due to multiuser interference similar to ergodic capacity 
for single-cell systems. 

2) Network Throughput: Define now the network throughput ratio gap QT to be the ratio of the 
throughput achieved by zero-forcing with perfect CSI to the throughput of finite rate feedback zero- 
forcing, which is given by QT = (1 + 135)^^^. The number of feedback bits per user for throughput ratio 
offset QT = r needs to scale according to 

Bqt > (M - 1) log2 /3 - (M - 1) log2 (rT^ -1^ bits/user, 

for any r satisfying 1 < r < (1 + (3)^^^. 
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Similarly to the transmission capacity offset, the number of feedback bits has to increase at least linearly 
with the number of antennas/streams and logarithmically with the target SINR constraint to maintain 
constant throughput loss. For instance, for /3 = 3 dB and M = = 4, at least B = 9 bits are required for 
a 3-dB offset. Furthermore, for r = 2^^\ the resulting feedback scaling takes on the following simple 

Af - 1 

Bqt = g /3dB bits/user. (24) 

Note that the feedback bit scaling takes on identical form with that in [19] for a 3-dB rate offset substituting 
transmit power P for target SINR /3. 



B. Feedback Rate Scaling 

In the previous section, we quantified the number of feedback bits required for a fixed gap degradation 
in the throughput and transmission capacity performance of zero-forcing with imperfect feedback. Here, 
we provide design guidelines on the feedback bit scaling for asymptotically vanishing performance loss 
due to partial CSIT. 

First, similar to MISO broadcast channels without inter-cell interference, both throughput and trans- 
mission capacity of zero-forcing with limited feedback are bounded with fixed B even if other system 
parameters grow large. Furthermore, if the feedback bits do not scale with M and/or (3, the throughput 
ratio and the transmission capacity offset become unbounded for asymptotically large values of M, /3. 

In the high antenna/stream regime, it can be shown that if the feedback load B is scaled with M at a 
rate strictly greater than (Af — 1) logg Af, i.e., B = (Af — 1) log2(Af'') for any r] > I, the transmission 

capacity offset converges to zero, i.e. lim AC = lim log(l + = 0, and the throughput 

A/-s>oo A/->oo T. 

-''A/ Sin ax 

ratio gap converges to one, i.e. lim QT = lim (1 + f36)^ — > 1. 

A/^oo A/— >oo 

The rate of convergence to zero and one respectively depends on r] and is faster with r] increasing. 
The throughput ratio converges to one also in the case where B scales superlinearly with the number 
of antennas, i.e. B = Af. Based on the above results, we establish that, at asymptotically high Af and 
under the aforementioned bit scaling, the network throughput and the transmission capacity of finite rate 
feedback zero-forcing converges to the perfect CSI throughput. In contrast, if the feedback rate is not 
properly adapted, the transmission capacity offset scales as AC = 0(Af^^^/"). 

In the high reliability regime (large (3), if the feedback load B is scaled with /3 at a rate strictly greater 
than (Af — 1) logg (3, i.e. B = Klogg /3 for any k > M — 1, the throughput ratio gap converges to one, i.e. 
lim QT = lim (1 + (36)^-^ 1. 

/3— i>oo (3^00 

Under the same bit scaling, the transmission capacity offset vanishes (asymptotically in (3), i.e. lim AC = 

Af-5-oo 

0. Thus, at large (3 (high SINR regime) and under the aforementioned bit scaling, the throughput of the 
finite rate feedback ZFBF converges weakly to the perfect CSI throughput. 



C. Optimal Number of Streams 

In this section, we investigate the optimal number of streams to be used per cluster in order to maximize 
the capacity. These results also provide useful insights on the feasibility and the potential gains of multi- 
stream, multiuser beamforming and adaptive beam selection in wireless ad hoc networks with imperfect 
feedback. 

For that, we consider again that K < M streams can be sent. We define the per-user throughput as the 
normalized network throughput over the number of users, i.e. 7^ = -^T. Taking the partial derivative with 
respect to the number of streams, we can show that 7^ is a decreasing function with K. This confirms 
the intuitive argument that, from the user perspective, employing single-stream beamforming {K* = 1) 
maximizes the per-user throughput. 
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However, the optimal number of users to serve may alter if we consider the system (broadcast cluster) 
overall throughput. In this case, there is a tradeoff between spatial reuse and feedback quality, i.e. for 
certain values of K and other system parameters, the spatial multiplexing gain may compensate for the 
performance degradation incurred due to finite rate feedback. For the network throughput, the complicated 
form of the interference constant Xk precludes a simple, closed-form expression for the optimal number 
of streams. Specifically, we have the following result: 

Proposition 1: The number of users to be served per broadcast cluster that maximizes the network 
throughput is 

K* = min (max ( \uj] , 1) , M) 

where u is the solution for K of 



K + ^ + log(l + = 1. 



Proof: The result following by taking the derivative of (ITOl ) with respect to K and finding the optimal 
value of K based on Fermat's theorem for the stationary points. ■ 
Remark 3: Using the large K approximation ~ 7rr(l — 2/a)Ka , we have that K* = x", where 
X is the solution of the polynomial equation Cix" + C2X^ — 1 = 0, where Ci = S-C^ _j_ log(l + /35) 
and C2 = A7rr(l — 2/a)Cmax- By the Abel-Ruffini theorem, no general algebraic solution exists for 
a > 5, however for this particular form of the polynomial equation, K* can be found in closed-form for 
a E {2, 3, 4, 5, 6, 8}; solutions for a = 6 and a = 8 will be the squared root of the a = 3 and a = A 
solutions, respectively. 

In Fig. Ul we plot the network throughput as a function of the node density for different number of 
streams. As indicated by the above analysis, fully loaded SDMA (K = M) is detrimental for throughput 
performance, while multi- stream is superior only in sparse networks. An adaptive scheme that sets the 
number of streams based on the large K approximation solution is also plotted. Interestingly, the adaptive 
scheme results in the optimal multi-stream transmission scheme within almost all range of intensity values 
A. Furthermore, as expected, throughput increases when feedback quality is improved. We also plot the 
scheme in which the number of streams are adapted according to the large K approximate value, which 
performs satisfactorily for low densities A. 

Regarding the transmission capacity, taking the derivative of (fT2l) with respect to K and finding the 
stationary points, we have 

Proposition 2: The optimal number of streams that maximizes the multi-stream transmission capacity 

is 

K* =min (max([z/] , 1) , M) 



where v is the solution of 



&X 



with C3 = + log(l + and C^ = - log(l - e) + log(l + 

Remark 4: Using the large K approximation, the optimal number of streams is given by 

r (l-2/a)CJ 
* = (2-2/„)C3 ■ 

In Fig. [8l we plot the multi-stream transmission capacity vs. outage constraint for different number 
of streams employed. For most relevant parameters, single-stream transmission is optimal, even for a 
large number of bits. Similarly to the throughput case, the adaptive scheme based on the simple large K 
approximate solution performs satisfactorily in a wide range of e values. 
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V. Conclusions 

We investigated the performance of zero-forcing precoding with limited feedback in single-hop ad hoc 
networks under a broad set of metrics and scenarios. The main takeaway of this paper is that SDMA may 
not be a wise use of transmit antennas in decentralized networks with both self and other user interference, 
as single-stream transmission maximizes both outage-based and average throughput in most practically 
relevant scenarios. In other words, a high density of single-stream communication links may be preferable 
than SDMA transmission with quantized channel state information. Our analytical framework enables 
us to quantify the effect of the residual multiuser interference due to quantized CSIT on network-wide 
performance and to properly adjust the number of streams and the feedback bit scaling to achieve certain 
level of rate performance. A key finding that the per-user feedback load must be increased almost linearly 
with the number of antennas and logarithmically with the target SINR. The techniques developed here are 
also relevant for the analysis of partial CSIT-based linear precoding in emerging heterogeneous network 
paradigms, including femtocells, relays, picocells, and WiFi hotspots. 

Further extensions to this work could include different limited feedback precoding schemes (e.g. MMSE) 
or how to exploit multiple receive antennas for interference cancelation. Future work could consider 
multi-hop networks with opportunistic routing and investigate potential SDMA gains into end-to-end 
performance (e.g. progress-rate-density). It would be also of interest to explore how opportunistic user 
selection affects the spatial reuse and rate, as well as how to properly design CSIT in networks with 
spatial randomness. 



Appendix 

A. Proof of Theorem [7] 

Let Cp denote the Laplace transform of the normalized aggregate interference from the Poisson field of 
interferers (inter-cluster) Ip = ^^^^(^x) Sik\Xi\~°' with fading marks Sik, defined as Cp{s) = E7 



e '^^p 



/o°° 6 ^^fj^{p)dp. The Laplace transform of the normalized intra-cluster interference due to multiuser 
transmission with quantized CSIT is denoted by Cq{s). Define also the random variable F = + J^. 



The outage probability is given by J'^''\l3k, a) 

pSokdk"' 



^('=)(/3fc,a) = l-P 



Ip + Ig + (t2 



>/3fc 



P {SINRfc > Pk}, which can be rewritten as 

1 - P {Sok > Mt(lp + 7, + aVp)} . 



The channel gain is given by Sok = |hfcWfc| 



\u ii2 it: |2 



|hfcf i3(l,M- 1), where B{1,M-1] 



is a Beta distributed random variable (r.v.) with shape parameters (1,M — 1) and independent of ||hfc| 
[|24l . The term \\hk\f is distributed as a chi-squared r.v. with 2M degrees of freedom denoted as xIa 
Thus, Sok ~ exp(l) (exponentially distributed with unit mean). Denoting (k 



f3kd'^, we have 



J'^'^Pk, a) - E [exp(-a(/p + /, + aVp))] 1 



(fe) 



_Cfe£f 

£p(Cfc)A(Cfe)e " 



(26) 



where step (a) is reached by conditioning on the aggregate interference Ip+Iq and (b) by the independence 
of the interference terms. 

The interferer marks in Ip are chi-squared distributed with degrees of freedom S^k = ||hjfcWj||^ ~ X2K 
since it is the sum of K i.i.d. exponential random variables. Thus, the Laplace transform for a Poisson 
shot noise process in with i.i.d. xIk distributed marks is given by [|27l 



Cp{s) 



E<5 



ie«(A) 



Si\Xi\-^) 



exp 



-A 



1 -E 



s 



'sS\x\ 



dx 



-As a Ik 



where 



2vr fK 



'-K 



E 



a ^ — ' \m 

m,=0 



B m + -,K 



m 



a 



a 



(27) 



(28) 
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with B{a, b) = f-^l - tf-^dt = ^J^g^ being the Beta function. 

For the marks in the interference term we have that |hfcWj |^ = ||hfc||^ |hfeWj|^. Denoting 0^ the angle 
between and h^, and decomposing the normalized channel vector as = (cos0fc)hfc + (smcpkWk, we 

2 ||2| |2 9 

have that |hfeWj| = ||hfc|| |vfeWj| sin 0^. Since wj is isotropic within the hyperplane and independent 
of Vfc, the quantity |vfeWj|^ is beta (1, M — 2) distributed, i.e. |vfcVi^j|^ ~ B{1, M — 2) [19] and independent 
of the quantization error sin0fc. Thus, the normalized intra-cluster interference can be rewritten as 

7, = d^°||h,f sin2 0, 13{1,M~2), (29) 

where ||hfe|^, 0^, and B{1, M — 2) are all independent. For the quantization cell approximation, the term 
X = ||hfc|| sin^0fc is a gamma distributed r. v. with shape M — 1 and scale S, i.e. X ~ Gamma(M — 1, 5). 
Therefore Y = XB{1,M - 2) is exponentially distributed with rate 1/6 BH. As 7g = d'^'^Z, where Z 
is the sum of {K — 1) i.i.d exponentially distributed r.v., i.e. Z ~ Gamma(ii' — 1,5), the normalized 
aggregate interference becomes Iq ~ Gamma(-ft' — 1, Sd'^°'). The Laplace transform of Ig is given by 

Substituting ^} and dSO]) in we obtain the resuk. 

B. Proof of Theorem 2] 

For a random variable X with probability density function fx {x) and cumulative distribution function 
(cdf) Fx{x), we have 

poo poo POO 1 rp / \ 

E{log(l+X)}= / log(l+x)/x(x)da:= / log(l + x)d [1 - Fx(a:)] — ^iidx 
Jo Jo Jo ^ + ^ 

where step (a) follows from integration by parts. The result is obtained by 

n, = / —-^dx = / j^—^f^-Y^.dx, with C, = and C, = AX^rf^ 

Jo 1 + Jo [l + x){l + dxy^^^ p 

C. Proof of Theorem |?] 

First, we consider the following upper bounds to the average rate achieved by the A;-th user: 

Tlu = E nog 1 + 4^' % " ^ U < E hog 1 ' ' ' " ^ 



h + h + ^Vp I ] { V ip + h + ^Vp 

= E jlog (^7p + 7, + ^ + IhfcWfcl' rf^"^ I - E jlog {l\ + ~h + -p 

< log {e |7^ + 7, + ^ + |hfcw,|' - E {log (7,) } 

< log(^E|7^}+E{74 + ^ + rf-"E{|hfcWfc|'}^ -E{log(d^"|hfcW,f)} (31) 



where Ip = Ip/p and Iq = Iq/ p are given in (|3]) and (HJ, respectively. In (a) we consider a lower bound 
to the inter-cluster interference, denoted as J^, and in (b) we apply Jensen's inequality to the minuend 
and neglect the noise and inter-cluster interference terms in the subtrahend. Step (c) follows from the 
independence of the interference terms and the received signal and by considering only one of the intra- 
cluster interference terms. 
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The normalized interference Ig ~ Gamma(-R' — 1, 5dj^"), thus we have E {lg\ = dj^"5{M — 1). The 

I 1 2 r I 1 2 ~i 

channel gain Ih^w^l is exponentially distributed with mean one, i.e. E jlh^w^l | = 1. In order to derive 

a lower bound, we neglect the contribution of the closest interferer to the inter-cluster interference. Using 

techniques from [10], the lower bound is given by = (vrArf^)^ ([|] — l) ^ . 

For the second term in (ISTT ). we need to compute 



E{log (|h,w,f )} = Ejlog (llhfcf |h,w/)} = E{log (||hfcf ) } + E {log (jh^w/)} . (32) 

For the channel norm we have ||hjt||^ ~ xIm^ thus Ejlog (||hfc||^)} = J^^logx ■ ^^^^^^y^dx = ^{M). 
Then, E jlog (^jhfcWjf ) | = E (log (sin^ 0^) } + E {log M - 2))}, with 

E{log(sinV.)} = log(x)dF,i„.^^(x) = 2^(M-l) log(x)x^^'^dx = - _ ^ 

and 

E{log(S(l,M-2))} = logx--^^^-^^^-^dx = -i7M-2 

where the cdf of the quantization error sin^ 0^ is given in [24J. By substituting the above quantities to 
(|3TI) and after some manipulations, we obtain (fTSi) . 



D. Proof of Lemma [7] 

The proof of this lemma is based on Holder's inequality [|26l . which states that for two measurable 
functions / and g defined on a Hilbert space S and 1 < p, q < oo with 1/p + 1/q = 1 

^Jfix)gix)\dx= (^Jjf ix)fdx^ (^j\g{x)\'^dx^ '\ (33) 

The result is obtained by applying twice (|33l) for the following decreasing real-valued and bounded 
functions: first /(x) = e-Cix^-C2xa ^^^-^ _ ]^ and then for f{x) = e~^^^ and g{x) = 

^-C'2xa g^^j^ gj^gp^ optimize p, q in order to obtain the tightest possible upper bound. Using the 
following forms 

Xi = / e-'^Mx = -, X2= / e-^^'^''dx = r"/2r(l + a/2), 



X, 



a 

^ dx n cscig7i){5 - l)i-^-55^S(c, ^ - 1) +2 Fi(l, c, 2 - i) 



{l + xy{l + 6x)9 Sg 

where csc(-) denotes the cosecant and 2Fi{x,y, z,w) the Gauss hypergeometric function, and after some 
algebraic manipulations, we obtain (fTSi) . 
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Fig. 1. Outage probability vs. node density at SNR = 20 dB, for a = 4, d = 1.5, and /? = 1 dB. 



SISO (M = 1) 

M = 2 (partial CSIT) 

M = 4 (perfect CSIT) 

M = 4 (partial CSIT) 



_ B= 12 



B=4 



0.09 
0.08 
0.07 
0.06 

I- 

t 0.05 

i 0.04 
1- 

0.03 
0.02 
0.01 


0.02 0.04 0.06 0.08 0.1 0.12 0.14 

Intensity, X 

Fig. 2. Throughput vs. node density for a = 4.2, d = 1.5, /3 = 3 dB, and SNR = 15 dB. 
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Fig. 3. Transmission capacity vs. outage constraint for a = 4.5, /3 = 1 dB, and SNR = 20 dB. 
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Fig. 6. Average rate vs. SNR for a = 4.2, d = 1, A = 0.05, M = 3, and B = 10 bits. 
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Fig. 7. Network throughput vs. node density for a = 4, d = 1.5, M = 4, B = 10 bits, /3 = 1 dB, and SNR = 15 dB. 
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Fig. 8. Multi-stream transmission capacity vs. outage constraint for a = 4.5, M = 4, B = 12 bits, P = 1 dB, and SNR = 20 dB. 



